2D to 3D transition in soap films demonstrated by microrheology 
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We follow the diffusive motion of colloidal particles of diameter d in soap films of varying thickness 
h with fluorescence microscopy. Diffusion constants are obtained both from one- and two-particle 
microrheological measurements of particle motion in these films. These diffusion constants are re- 
lated to the surface viscosity of the interfaces comprising the soap films, by means of the Trapeznikov 
approximation [A. A. Trapeznikov, PICSA (1957)] and Saffman's equation for diffusion in a 2D fluid. 
Unphysical values of the surface viscosity are found for thick soap fllms {h/d > 7 ± 3), indicating a 
transition from 2D to 3D behavior. 



PACS numbers: 47.57.Bc, 68.15.+e, 87.16. D-, 87.85.gf 

A soap film is a thin layer of fluid stabfized by two 
surfactant layers that buffer it from air phases above and 
below. In the early 18th century, Sir Isaac Newton mea- 
sured the thickness of the fluid layer to ^ 10 nm precision 
[l|. In fact, because of the similarity of thin "Newton 
black fllms" to planar lipid bilayers, soap fllms have been 
proposed as models for cell membranes [2] . The analogy 
with membranes extends to considering a thin soap fllm 
as a 2D fluid [3|]. This has motivated the use of soap 
fllms to study turbulence in 2D [J, Q , as well as inform- 
ing the physics of drainage in foams [g|. However, soap 
fllms have a nonzero thickness, and presumably under 
some conditions the model of the fllm as a 2D fluid is 
inappropriate. 

A previous study by Cheung et al. [3|] quantifled the 
hydrodynamics of a single soap fllm for the special case 
where the diameter d of embedded tracer particles was 
assumed to be the same as the thickness h of the film. 
They observed the relative diffusive motion of pairs of 
these particles, and found that this relative diffusion de- 
pended logarithmically on the separation between the 
particles, indicating 2D fluid-like behavior. Clearly, for 
thicker films where h ':$> d, 3D behavior must be recov- 
ered. This transition from 2D to 3D behavior has not 
been demonstrated in any study to date. 

In this Letter, we use the thermal motion of embed- 
ded polystyrene particles to study soap fllms of varying 
thickness h, to clarify which aspects of particle motion 
arise from 2D hydrodynamics and which from 3D hydro- 
dynamics. For small particles in thick films {h/d > 7), 
the measured particle diffusivity is similar to free 3D dif- 
fusion in the fluid comprising the film. For thin films 
{h/d < 7), particles diffuse noticeably faster, suggesting 
that particle drag is more due to 2D hydrodynamics with 
an effective 2D viscosity. Measurements of the correlated 
motion of pairs of particles show that all soap films have 
2D-like long-raiige correlations. The classic Trapeznikov 
approximation Q connects the 2D and 3D properties of 
the film by modeling the soap fllm as a 2D interface with 
an effective surface viscosity 77s ^eff given by 

r]s,cS = r]hu\kh + 2r]int (1) 

in terms of the 3D viscosity 77buik of the fluid in the film. 



and the 2D surface viscosity rji^t of the surfactant layers; 
see Fig. 1. Our results show that Eq. (1) and Fig. 1 
are valid for thin fllms but not for thick fllms where 3D 
hydrodynamics becomes important. These observations 
lead us to conclude that a transition from 2D fluid to 3D 
bulk behavior occurs at around h/d « 7 ± 3, the first 
experimental demonstration of such a transition. 
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FIG. 1: Schematic of the Trapeznikov approximation, where 
the entire soap film is approximated as a single interface in 
contact with bulk air phases. 

We use mixtures of water, glycerol and the commer- 
cially available dishwashing detergent Dawn to prepare 
our soap films. This particular brand was chosen as it 
has proved quite successful at making long-lasting soap 
fllms Q; further, the preponderance of its use in the lit- 
erature allows us to compare our results with previous 
work [3, [y, Q . The concentration of Dawn is kept fixed 
at 2% by weight in our soap solutions to try to maintain 
a constant interfacial viscosity rj-mt for all fllms. The fluid 
viscosity r/buik in our films is controlled by changing the 
ratio of water and glycerol in the soap solutions. Fluo- 
rescent polystyrene spheres (Molecular Probes, carboxyl 
modified, d — 210 or 500 nm) are added to these mixtures 
to act as tracer particles. 

Stable soap fllms are created by dipping and drawing 
out a circular stainless steel frame of thickness 1 mm 
from the soap solutions. The steel frame is enclosed in 
a chamber that minimizes convective effects in the soap 
fllm while maintaining its relative humidity. We then 
image the particles by fluorescent microscopy; soap films 
containing 500 nm particles are imaged with a 20x ob- 
jective (numerical aperture — 0.4, 465 nm/pixel) while 
those with 210 nm particles are imaged with a 40x ob- 
jective (numerical aperture — 0.55, 233 nm/pixel). For 
each sample, short movies of duration ~30 s are recorded 
with a CCD camera that has a 640 x 486 pixel resolu- 



tion, at a frame rate of 30 Hz. After each movie, we 
transfer the film to a spectrophotometer and its thick- 
ness h is determined from the transmitted intensity [lO] . 
The movies are then analyzed by particle tracking to ob- 
tain the positions of the tracers. From the particle po- 
sitions, we determine their displacements by the relation 
Ar(i, t) = r{t+T)—'r{t), where t is the absolute time and 
r is the lag time. Finally, any global motion is subtracted 
from these displacements to eliminate the remnant effects 
of convective drift caused by the air phases that contact 
the soap film. 
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FIG. 2: The solid symbols are the mean square displacements 
for five soap films with increasing h/d ratios. See Table I for 
the viscosities r^buik that correspond to these films. The open 
symbols represent data from the bulk (3D) solutions used to 
make the soap films for four cases (b, g, i/j). The dashed line 
with a slope of 1 is a guide to the eye. 

To quantify particle motion, we use the vector displace- 
ments of the particles to calculate the mean square dis- 
placements (MSDs) (Ar^) as a function of the lag time t. 
Figure 2 shows measurements of (Ar^) for five different 
soap films, where h/d ranges from 0.6 to 15. The viscos- 
ity r^buik of the fluid layer of these soap films is given in 
Table I. At long lag times, all the MSDs are linear with 
respect to r, irrespective of h/d and rybuik- This is ex- 
pected, as both the fluid layer and the interfaces compris- 
ing the soap fllm are viscous. We extract a one-particle 
self-diffusion coefficient, I?s,ip, from the measurements 
according to the equation (Ar^) = ArDg^ip- Further, 
from Table I and Fig. 2 it is clear that increasing both 
??buik and h/d tends to slow the diffusion of the parti- 
cles. Films made from more viscous bulk fluids tend 
to be thicker (see Table I), and so these single-particle 
measurements do not clearly distinguish the influence of 
'7buik and h/d, although it is obvious that increasing yybuik 
should slow diffusion, and plausible that increasing h/d 
might also slow diffusion. This latter effect is suggested 
by comparing films i and j in Fig. 2, which have the same 
^ybuik! the motion is slower for the thicker film j (upward 
solid triangles). A further suggestion that thicker films 



TABLE I: Parameters for all the soap films described in this 
paper. 77buik (determined from diffusivity measurements in 
bulk solutions) has an error of ±5%, and values of h and d 
are certain to within ±2%. The uncertainities in rjint, derived 
from Eq. 1 and 2, are given in the brackets. 
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h 


d 


77int(lp) 


'?int(2p) 


(mPa-s) 


(nm) 


(nm) 


(nPa-m-s) 


(nPa-m-s) 


a. 2.0 [3] 


400 


400 


0.20 (± 0.03) 


0.47 (± 0.06) 


b. 2.3 


305 


500 


0.63 (± 0.06) 


1.02 (± 0.10) 


c. 3.0 


640 


500 


0.49 (± 0.09) 


0.62 (± 0.12) 


d. 6.0 


510 


500 


0.89 (± 0.2) 


0.84 (± 0.2) 


e. 10.0 


1340 


500 


0.34 (± 0.5) 


2.26 (± 0.7) 


f. 25.0 


1100 


500 


-0.30 (± 0.9) 


4.35 (± 1.5) 


g. 10.0 


780 


210 


0.12 (± 0.26) 


1.64 (± 0.4) 


h. 25.0 


2184 


210 


-8.92 (± 1.3) 


27.2 (± 4) 


i. 30.0 


2100 


210 


-10.6 (± 1.5) 


25.0 (± 4) 


j. 30.0 


3000 


210 


-15.5 (± 2.1) 


65.2 (± 8) 



result in slower diffusion comes from comparing the mo- 
tion within the films (solid symbols in Fig. 2) with motion 
in the 3D fluid solutions the films are made from (open 
symbols in Fig. 2). For a thin film {h/d ~ 0.6, stars) the 
particle motion is much faster in the soap films than in 
the corresponding 3D solution. For the thickest films we 
study {h/d w 10 — 15, solid diamonds and triangles) the 
motion in the soap film is comparable to the motion in 
the corresponding 3D solution (open diamonds). 

To further understand the hydrodynamics and how 
particle motion compares in thick and thin filrns, we 
use the correlated motions of particles [ll|, [I^, \lM to 
probe flow fields in these soap films. Bricfiy, we look 
at the product of particle displacements Drr{R,T) — 
{/\r\.{t,T)Arl{t,T)5{R — R^^ {t)))i^j^t, where i, j are par- 
ticle indices, the subscripts r represent motion parallel 
to the line joining the centers of particles, and R^^ is 
the separation between particles i and j. Similar to [l3|, 
we observe D^r ~ t, which enables the estimation of a 
T-independent quantity (J)„/t)ti depending only on R 
and having units of a diffusion constant. 

In Fig. 3 we show (J)„/t) as a function of the separa- 
tion R for the five soap films described in Fig. 2, with an 
additional data set from Ref. Q included. The motion 
of a tracer particle creates a flow field in the soap film 
that affects the motion of other particles, and the cor- 
relation function indicates the spatial extent of this flow 
field. The dashed line in Fig. 3 represents the form of 
the correlation function in a 3D fluid [ll| ; therefore it is 
clear that the motion is correlated over larger distances 
in soap films than in 3D. This long-ranged behavior is 
characteristic of 2D fliuds [IJ, [l^. Further, similar to 
the trend seen in the MSDs, increasing h/d and T^buik 
lowers the value of the correlation function {Drr/T) for 
the same separation R. As {Drr/T) dimensionally repre- 
sents a diffusion constant, slower diffusion decreases its 
magnitude. Finally, the correlation functions for all six 
soap films are similar in shape. This is evident from the 
form of the function, Ahi R + B that has been used to 
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FIG. 3: Two-particle longitudinal correlation function 
(Drr/T) for the five soap films described in Fig. 2, symbols 
being the same. An additional data set from ^ has been in- 
cluded (open hexagons). Solid lines are empirical fits to the 
data of the form A In R+B. The dashed line is the form of the 
correlation function expected for a 3D fluid with Jjbuik = 6.5 
mPa-s (squares). 
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empirically fit all the six curves. 

The data set from Cheung tt al. Q (open hexagons) 
requires explanation. In their paper, the data was pre- 
sented as a two-particle MSD, (Ai?^) = {{[rj{t + r) - 
ri{t + t)] — [rj{t) — ri{t)]}'^), which measures the rel- 
ative diffusion between particles i and j. This was 
done for a fixed r (=1/30 s), and the resulting rela- 
tive diffusion decomposed into two components (Ai?^) ~ 
(Ai?|) -I- {AR^}, representing displacements parallel and 
perpendicular to the lines joining the centers of the par- 
ticles. It is then straightforward to show that (Drr/T) = 
i2{Arf) - (Ai?^))/2r « ((Ar^) - (Ai?f))/2r (see M)- 
We plot the data in this form in Fig. 3, evaluated at r = 
1/30 s. 

The data shown in Fig. 3 can be used to extract 
a single-particle self diffusion constant, Ds2p, but now 
measured from two-particle correlations. This is done 
by extrapolating the correlation functions (Drr/T) to 
R = d/2. The reason for this choice is that the single 
particle diffusion constant must be recovered from the 
two-particle measurement when extrapolated to the par- 
ticle radius [111]. We then deduce that Drr{R — c?/2, r) = 
(Ar^) ss 2Ds^2pT [13 and use the fitting functions shown 
in Fig. 3 to determine the value of Ds,2p for each soap 
film. However, this extrapolation process has limitations: 
nearby particles at interfaces can have strong interac- 
tions, either electrostatic or through capillary forces. Our 
particle concentration was chosen to have spheres be no 
closer than R ~ 5d, to avoid such effects. The excellent 
agreement between the fitting function and all our data 
gives us confidence, however, that the extrapolation is 
valid even when R < 5d. 

Insight into the transition from 2D to 3D behavior can 



FIG. 4: (a) Two-particle diffusion constant Ds,2p plotted 
against the one-particle Ds,ip for the six soap fllms described 
in Figs. 2 and 3. Four additional data sets have been in- 
cluded, details of which are given in Table I. The straight line 
indicates equality between the diffusion constants. Errors in 
Ds,ip and Ds,2p axe ±5%, similar to that in 77buik, and are 
smaller than the size of the symbols in the flgure. (b) Interfa- 
cial viscosity j^int as a function of h/d. Inset: Magnified view 
of 7?int for h/d < 4, with error bars included. 



be gained by plotting the diffusion constants obtained 
from the two methods against each other, as shown in 
Fig. 4(a). The two diffusion constants agree with each 
other for thin soap films while for thicker films, the devia- 
tion from equality lies beyond experimental error (points 
h, i, and j). These results can be interpreted using the 
Trapeznikov approximation described in Fig. 1, which 
models the soap film as a 2D interface with an effective 
surface viscosity ?7int given by Eq. (1). Based on this 
approximation, our system reduces to that of a particle 
diffusing at an interface in contact with bulk air phases. 
The diffusion of a disk or sphere IitII at such an interface 
has been described by Saffman [l8| and is given by 
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where Tyair is the viscosity of the bulk air phases and ^e 
is Euler's constant. Equation (2) holds if 27ys.off S> '7airi^, 
which is true for our data due to the low value of the air 
viscosity (r^air = 0.017 mPa-s). 

We now attempt to determine the interfacial viscos- 
ity of the surfactant layers by using Saffman's equation 
[Eq. (2)] to convert measurements of -Ds.ip and Ds^2p 



into rjs^oS, and then using the Trapeznikov approxima- 
tion [Eq. (1)] to determine rji^t = l/2(r?s,eff-??buik/j)- For 
all films, we average Dg^ip and Ds,2p to determine ryint, 
which is plotted in Fig. 4(b) as a function of h/d. For the 
thin films the interfacial viscosity shows roughly constant 
behavior while for thick films the variation in T^int is quite 
pronounced beyond the experimental uncertainty in the 
measurements (see Table I). The inset to Fig. 4 (b) shows 
a magnified view of the interfacial viscosity for h/d < 4, 
where it is clear that 77int is nearly constant with an aver- 
age value of 0.97±0.55 nPa-s-m. This is expected because 
the same concentration of Dawn surfactant has been used 
in all these soap films. 

For thick films, the one-particle measurements Ds^ip 
give large negative values of r/int (refer Table I) imply- 
ing that the single particle diffusivities are significantly 
faster than that predicted by Eqs. (1,2). From Fig. 2, it is 
clear that the 3D Stokes-Einstein equation for diffusion, 
Ds,ip = ksT /{Ziirjx^-aw^d), is sufficient to explain the mo- 
tion of the probe particles in thicker soap films, without 
the need to invoke Saffman's equation. This makes sense, 
as in the limit of a 3D system {h -^ cxd), Eq. (1) predicts 
an infinite surface viscosity, which has no physical mean- 
ing. The apparent negative values of rj-mt for h/d > 7 
indicate that the 3D limit is already evident for films of 
this thickness. In contrast with the one-particle measure- 
ments, the two-particle measurements in thick films give 
large positive values of yyint (see Table I), again contra- 
dictory to the low values determined in thin films. An 
alternate way to state this is that the Trapeznikov ap- 
proximation predicts an effective surface viscosity that is 
too small, if we use r^jnt based on the thin film measure- 
ments. 

This leaves us with a puzzle; even for these thick films 
the two-particle correlation functions are long-ranged in- 
dicating that the soap films behave like a 2D fluid. In 
fact, the behavior of the correlation functions as a func- 



tion of R for all soap films can be explained by consid- 
ering the following. Locally, the films likely behave as a 
3D fluid [l^l- We hypothesize that the correlation func- 
tions in the thick films would then decay as 1/R at very 
short separations (d/2 < R < h) but more slowly at 
larger separations {R > h) . The extrapolation of Drr to 
R = d/2 thus underpredicts -Ds.2p, explaining the over- 
estimation of 7/int for the thick films. At intermediate 
separations, because of conservation of fluid momentum, 
all soap films behave as 2D fiuids [I^. Therefore, the 
correlation functions decay in a logarithmic fashion for 
those separations, as seen by the form of the fitting func- 
tions in Fig. 3. However, the logarithmic divergence of 
the correlation function is cut off at a length scale where 
stresses in the air phase from motion of the the tracers be- 
come important. This length scale, related to ?7s,cff/'7air, 
is the separation at which the correlation functions begin 
to decay more rapidly, indicating a final crossover to 3D 
fiuid like behavior. 

Our work describes a transition from 2D to 3D behav- 
ior when the thickness of soap films are changed with 
respect to particle size, at a ratio of h/d « 7 ± 3. This 
demonstrates that the hydrodynamics of thin films de- 
pend on the ratio of the dimensions of the film and the 
probe particle, rather than intrinsically on the film it- 
self. The particular value of h/d « 7 ± 3 is an empirical 
determination of when 3D shear gradients in the fiuid 
layer dominate dissipation of stress from the motion of 
the probe, in comparison to the air phase. For thin films 
showing 2D behavior, the air phase is crucial for stress 
dissipation, as demonstrated by the presence of ?7air in the 
Saffman equation [Eq. (2)]. In all our films, the interface 
is relatively mobile, that is, ?7int is small when compared 
to the contribution from the fluid layer yybuik^- Changing 
the nature of the interface, such as making the interface 
more rigid {rjint ^ ''7buik^, for instance) will likely change 
where the transition occurs. 
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